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ABSTRACT 

We present the results of an XMM-Newton observation of the radio-quiet X-ray pulsar IE 1207-5209 located 
at the center of the shell-like supernova remnant G 296.5H-10.0 . The X-ray spectrum is characterized by the 
presence of two phase-dependent absorption lines at energies ^^0.7 keV and ^1.4 keV. Moreover, these broad 
spectral features have significant substructure, suggesting that they are due to the blending of several narrower 
hnes. We interpret such features as evidence for an atmosphere containing metals and a magnetic field value of a 
few 10^^ G, consistent with the observed spin-down rate P=(1.98±0.83)x lO"'"* s s"'. Since IE 1207-5209 is the 
only X-ray emitting pulsar showing evidence of such features, we tentatively link them to the unique combination 
of age and energetics that characterize this object. We suggest that a young age and a low level of magnetospheric 
activity are favorable conditions for the detection of atomic spectral features from Z > 1 elements in neutron star 
atmospheres, which would be either blanketed by a thin layer of accreted hydrogen in older objects or masked by 
non-thermal processes in young energetic pulsars. 

Subject headings: Stars: neutron; Pulsars: individual (IE 1207-5209 ); X-ray: stars; 



1. INTRODUCTION 

The thermal emission from the surface of a neutron star 
traces the star's cooling history. Its study can thus provide in- 
valuable information on the poorly known equation of state of 
matter at super-nuclear densities and on physical processes in 
strong magnetic fields. Satellite observations carried out in the 
last decade have clearly shown the thermal origin of the soft 
X-ray emission ('^0.1-3 keV) from a handful of middle-aged 
radio pulsars (Becker & Triimper 1997) as well as from several 
radio-quiet neutron stars (Caraveo, Bignami & Triimper 1996; 
Treves et al. 2000). 

It is expected that the presence of an atmosphere on the neu- 
tron star surface distort the emerging radiation by altering the 
blackbody energy distribution and introducing absorption fea- 
tures (see, e.g., Zavlin & Pavlov 2002 for a recent review). Fits 
with atmospheric models generally yield emitting regions com- 
patible with standard neutron star dimensions, thus providing 
indirect evidence for the presence of an atmosphere. However, 
until recently no convincing evidence for the absorption fea- 
tures predicted by these models was found, thus leaving sub- 
stantial uncertainty on the atmospheric composition and mag- 
netic field. 

Here we report on an XMM-Newton observation of the radio- 
quiet neutron star IE 1207-5209 . Previous X-ray, optical 
and radio observations (Bignami, Caraveo & Mereghetti 1992; 
Mereghetti, Bignami & Caraveo 1996; Vasisht et al. 1997) 
strongly suggested a neutron star nature for this source, located 



close to the geometrical center of the shell-like supernova rem- 
nant G 296.5H-10.0 (Roger et al. 1988). This was confirmed by 
the discovery of fast X-ray pulsations with period f =0.424 s 
(Zavlin et al. 2000). A subsequent measurement (at the ^ 2a 
level) of a positive period derivative P = (2.0^} 3) x 10"'"* s s"' 
(Pavlov et al. 2002a) results in a large discrepancy between the 
pulsar's characteristic age Tc = -^ = 200-900 kyrs and the age 
of 7 kyrs (with a factor 3 uncertainty) estimated for the associ- 
ated SNR (Roger et al. 1988). 

Our XMM-Newton data show the presence of broad absorp- 
tion features at ~0.7 and ^1.4 keV in the spectrum of IE 1207- 
5209 . Such features have been independently discovered in 
data from the Chandra satellite (Sanwal et al. 2002). The high 
throughput of the XMM-Newton telescope, coupled to the good 
spectral and timing resolution of the European Photon Imaging 
Camera (EPIC) instrument, allow us to show that the lines have 
a significant substructure, which varies with the phase of the 
pulsar. 

2. DATA ANALYSIS AND RESULTS 

The XMM-Newton observation of IE 1207-5209 started on 
December 23, 2001 at 19:13 UT and lasted 28.4 ks. We con- 
centrate here on data obtained with the EPIC instrument, which 
consists of two MOS CCD detectors (Turner et al. 2001) and a 
PN CCD instrument (Striider et al. 2001), for a total collecting 
area > 2500 cm^ at 1.5 keV. 

The PN camera was operated in Small Window mode in or- 
der to have a time resolution (6 ms) adequate to study the pul- 
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2.2. Spectral Analysis 

All spectral modeling was done with XSPEC VI 1 and using 
the most recent EPIC response matrices ' . For the FN analysis 
we used both single and double events in the 0.3-3 keV energy 
range (we checked a posteriori that the results do not change 
by using only single events). The spectra were binned to have 
at least 40 counts per channel and to over-sample by a factor 3 
the instrumental energy resolution. 

As found by Sanwal et al. (2002), fits with single-component 
models (power law, thermal bremsstrahlung, blackbody) gave 
unacceptable results (x^/dof > 4.4, dof=81) due to the pres- 
ence of broad absorption features at ~0.7 and ~1.4 keV. For 
completeness, we also tried the atmospheric models available 
in the XSPEC spectral fitting package, although they refer to 
weakly magnetized neutron stars (B < 10^- lO'" G) and might 
not be appropriate in the case of IE 1207-5209 , where the 
pulsations and the measured P testify the presence of a higher 
magnetic field (B~ 3 x 10'^ G). Hydrogen-atmosphere models 
(Zavlin, Pavlov & Shibanov 1996), which do not predict fines in 
the observed energy range, gave residuals very similar to those 
of the blackbody model, but, as expected, a lower temperature 
(kT~0.09 keV, as measured by an observer at infinity) and a 
larger emitting region. Models implying atmospheres with so- 
lar abundance or iron composition (Gansicke, Braje & Romani 
2002) gave unacceptable fits: the absorption lines they predict 
are too narrow to reproduce the features visible in the EPIC 
spectra. 

1E1207-5209 



sations without sacrificing the imaging, while both MOS CCDs 
were in the Full Frame mode (2.6 s resolution). All the de- 
tectors used the Medium thickness filter. AH the data were 
processed with the XMM-Newton Science Analysis Software 
(SAS Version 5.3). After screening to remove time intervals 
of high proton background which affected the MOS data, and 
correcting for the dead time, we obtained net exposure times 
of 18.7 ks, 22.3 ks and 24.8 ks in the PN, MOSl and MOS2 
respectively. 

IE 1207-5209 was clearly detected at a position consistent 
with previous measurements (Mereghetti et al. 1996), with 
net count rates (0.3-3 keV) of 1.352±0.008 counts s"^ and 
0.362±0.004 counts s~\ respectively in the PN and in each 
MOS. 

A circular extraction region with radius 40" was used for the 
timing and spectral analysis of the PN data. The background 
shows a slight intensity gradient across the chip. We verified 
that using background regions at different positions did not af- 
fect significantly the best fit parameters (the source count rate is 
~30 times greater than that of the background in the 0.3-3 keV 
range). We therefore finally used for the background spectral 
extraction a box of dimensions 4' x 2' located to the north of 
the source. 

For both MOS cameras we used a circular source extraction 
region (radius 40") and the background spectrum was estimated 
from a concentric annulus with radii 75" and 200". This re- 
gion is entirely contained in the central chip and is not affected 
by contamination from the internal Si-K fluorescence, largely 
present near the CCD edges. 

2.1. Timing Analysis 

For the timing analysis we used only the PN counts with pat- 
tern 0-4 and with energy in the range 0.2-2.5 keV. The times 
of arrival were converted to the Solar System Barycenter and 
folded, with 8 phase bins, in a range of trial periods around the 
expected value. This gave a significant detection of the pul- 
sation with a maximum x^=110 at P=424.131 ms. To deter- 
mine more accurately the period value we fitted the versus 
trial period curve with the appropriate sinc^ function, and com- 
puted the error on P using the relation between maximum 
and uncertainty derived by Leahy (1987). This resulted in our 
best estimate of P = 424.13084±0.00046 ms. For an indepen- 
dent assessment of the period uncertainty, we generated arti- 
ficial data sets with the same properties (duration, number of 
counts, pulsed fraction, etc..) of our observation and analyzed 
them in the same way. The difference between the derived and 
the true period was found to have a Gaussian distribution with 
a = 5 X 10"^ s, thus confirming the above estimate of the period 
uncertainty. 

Comparison with the period measured in January 2000 with 
Chandra (as reported in the re-analysis of Pavlov et al. 2002a) 
yields a period derivative P=(1.98±0.83)x 10"'"* s s"'. 

The folded fight curve is nearly sinusoidal, with a ~8% 
pulsed fraction (Fig. 1) and no significant energy dependence. 
We have also examined the profiles in the 0.3-1 keV and 1-1.7 
keV energy bands (16336 and 7639 counts, respectively) with- 
out finding the phase shift reported by Pavlov et al. (2002a). 
Simulations show that a phase shift of ~0.4-0.5 has a prob- 
ability much smaller than 1% to be obtained by chance with 
the available statistics. This suggest that the pulsar fight profile 
might be time variable. 

' PN version 6.1 epn_sw20_sdY9_me(lium, epn_sw20_sY9_medium; MOS ml_med_v9q20t5r6_all_15.rsp, m2_med_v9q20t5r6_all_15.rsp 




Phase 

Fig. 1 . — Light curve of IE 1207-5209 for three different energy ranges. 

Acceptable fits could only be obtained by adding to the mod- 
els two absorption lines, which for simplicity we modeled with 
Gaussian profiles. As shown in Table 1, the line parameters 
are only sUghtly dependent on the model adopted for the con- 
tinuum. The single blackbody yields a lower interstellar ab- 
sorption, while the other models give higher values, closer to 
the estimates obtained from radio observations (Giacani et al. 
2000). Although the blackbody plus power-law gives formally 
the best fit (see Fig. 2), we beUeve that this simply reflects 
the shortcomings of other models alone to reproduce the low 
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energy part of this complex atmospheric spectrum, rather than 
being evidence for a distinct non-thermal component. 

The fit residuals (lowest panel of Fig. 2) show that the broad 
lines have a profile characterized by the presence of significant 
sub-structure. Deviations from a smooth gaussian profile are 
seen at ~0.6 keV, ~0.75 keV and ~0.85 keV. Although the 
exact significance of such features is difficult to quantify, this 
may imply that the broad absorption below 1 keV is produced 
by blending of several narrower lines. 

Another interesting feature is shown by the residuals near 2 
keV. This was also noticed in the Chandra spectrum by San- 
wal et al. (2002), who could not exclude an instrumental effect. 
The fact that the same line is possibly present in the XMM- 
Newton data, even if at low significance, indicates that it might 
be really present in the source. 

We performed a similar spectral analysis based on the data 
from the two MOS cameras. This led to results entirely consis- 
tent with the ones discussed above, but with larger uncertainties 
caused by the lower counting rate. All the strongest spectral 
features were clearly visible in both CCDs. 

Blackbody + Power Law fit 
EPIC PN spectrum 



Structure, suggesting that several narrow lines contribute differ- 
ently at the various phases. 
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Fig. 2. — Fit with a blackbody plus power law and two gaussians. In the up- 
per panel the data are compared to the model folded through the instrumental 
response. The middle panel shows the residuals in units of sigma, the lowest 
one shows the residuals obtained by removing the lines from the model. 



2.3. Phase-resolved spectroscopy 

To search for possible phase-dependent spectral variations, 
we divided the PN data in four sets corresponding to the phase 
intervals indicated by the vertical lines in Fig. 1 . The result- 
ing spectra were fitted with blackbody models, keeping the 
absorption fixed at the best fit value of the average spectrum 
(N^f=3 X 10^" cm"^). While similar temperatures were obtained 
in the four phase intervals, the fit residuals (Fig. 3) clearly indi- 
cate that the absorption features are phase-dependent. In partic- 
ular the line at ~ 1 .4 ke V is virtually absent at the pulse peak and 
is more pronounced during the minimum and the rising parts of 
the pulse profile (see Table 2). Shape variations are also visi- 
ble for the lower energy feature, which shows a variable sub- 
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Fig. 3. — The four panels show the residuals of blackbody fits with N/f 
fixed at 3x 10^'^ cm~^ for the phase intervals shown in Fig. 1. 



3. DISCUSSION 

Our XMM-Newton observation confirms the presence of 
spectral features in the spectrum of IE 1207-5209 , as re- 
cently reported by Sanwal et al. (2002), and allows us to study 
them in more detail. In particular, both the PN and the two 
MOS spectra indicate that the feature at ~0.7 keV cannot be 
well described by a single gaussian line due to the presence of 
significant substructure. The PN spectra show phase-dependent 
variations, particularly pronounced for what concerns the inten- 
sity and width of the line at ^ 1 .4 keV. The presence of phase- 
dependent variations, leads to the obvious, but important con- 
clusion that the lines are linked to the rotating neutron star and 
not due to absorption in the line of sight. 

According to Sanwal et al., an interpretation of these spectral 
features in terms of cyclotron resonance lines is difficult, con- 
sidering the magnetic field inferred from the P measurement, 
B~3 X 10^^ G, and the relative intensity of the two lines. This 
view has been criticized by Xu et al. (2002), who considered 
the possibility that IE 1207-5209 be a bare strange star, spin- 
ning down in the propeller regime due to the presence of a fossil 
disk. 

Considering the conventional scenario of an isolated neutron 
star, a more likely possibility, also suggested by the structured 
and phase-dependent profiles seen with EPIC, is that these fea- 
tures result from atomic transitions. Pavlov et al. (2002b) in- 
terpret them as Hell lines in a strong magnetic field (B^(1.4- 
1.7)x 10'"* G), and derive for the neutron star a radius to mass 
ratio of (8.8-14.2) km Mq . However, such a magnetic field is 
much higher than the value (3±0.6) x 10^^ G estimated, assum- 
ing magnetic dipole braking, from the spin-down value (now 
confirmed by the XMM-Newton data). This implies either that 
the observed P is affected by ghtches and/or significant timing 
noise (and the true spin-down is of the order of P^ 6 x 10"^^ s 
s"') or that field components stronger than the dipole are present 
on (some region of) the neutron star surface (e.g. due to an off- 
centered dipole). 

Alternatively, the Unes can be due to heavier elements in 
a more conventional magnetic field. For example, models of 



4 



magnetized (B 10'^ G) iron atmospheres (Raj ag opal, Ro- 
mani & Miller 1997) predict, for temperatures ~ 10^ K, many 
absorption lines due to atomic transitions in the range above 
0.3 keV. Owing to the magnetic field and temperature varia- 
tion across the neutron star surface, these lines will probably 
be blurred. Indeed, the results of our phase-resolved spec- 
troscopy show that different physical conditions, most Ukely 
due to changing magnetic field configurations, are present on 
the neutron star regions responsible for the emission visible at 
different phases. As a consequence, the physical parameters in- 
ferred from the fit to phase-averaged spectra should be taken 
with some caution. It is also likely that the presence of several 
atomic absorption lines and edges affect the the parameters of 
the continuum as derived from medium resolution spectra. ^ 

IE 1207-5209 is the only neutron star in which significant 
X-ray absorption features have been detected so far (excluding 
of course the bright neutron stars accreting in binary systems). 
Observations with high statistics and good spectral resolution 
have been recently carried out for several thermally emitting 
neutron stars of various classes: the middle aged radio pulsar 
PSR B0656H-14 has a spectrum well described by two black- 
body components, without spectral lines in the 0.15-0.8 keV 
range (Marshall & Schultz 2002). No lines were found in the 
8.4 s radio-quiet pulsar RX J0720.4-3125 (Paerels et al. 2001) 
and in the nearby isolated neutron star RX Jl 856.5-3754. The 
latter was observed for more than 500 ks with Chandra and its 
spectrum was found to be well described by a blackbody func- 
tion with temperature kTBB=61 eV without evidence for any of 
the spectral fines or edges predicted by the models (Drake et al. 
2002). No fines have been detected in the spectra of the Vela 
pulsar (Pavlov et al. 2001a) and of the millisecond pulsar PSR 
10437^7 15 (Zavlin et al. 2002), which are comparable in qual- 
ity and statistics to the data presented here. Why is IE 1207- 
5209 different from the other coofing neutron stars which have 
been deeply scrutinized for the presence of lines with negative 
results? 

Before trying to answer this question, we must address the 
problem of the age of IE 1207-5209 . In this respect, its as- 
sociation with G 296.5+10.0 is crucial. This SNR has a re- 
markable bilateral symmetry. The two radio arcs that compose 
the shell have different curvatures, indicating a larger expansion 
velocity for the ejecta on the western side. The likely site of the 
supernova explosion is thus relatively well determined to lie to 
the east of the geometric center, within 8' from the position of 
the pulsar. The spatial coincidence between the pulsar and the 
SNR is remarkable, also in view of the relatively small number 
of such objects at this galactic latitude (b=10°). Furthermore, 
HI observations of this region (Giacani et al. 2000) indicate a 
spatial correlation between the two objects also for what con- 
cerns their distance. Thus we can consider this as one of the 
strongest neutron star/SNR association. An age of 7 kyrs has 
been derived for G 296.5+10.0 (Roger et al. 1988), very dif- 
ferent from the pulsar characteristic age Tc = (340 ± 140) kyrs. 
Although the SNR age estimate depends on several uncertain 
parameters, like its linear size, initial kinetic energy (£51 x 10^' 
ergs) and the density of the interstellar medium (n cm"^), it is 
very unlikely that G 296.5+10.0 is older than a few x lO"* years. 
For example, even taking the maximum distance of 4 kpc com- 
patible with the HI observations (Giacani et al. 2000), an age of 
^100 kyrs can be obtained only for (n/£5i) ^ 1. This would 

^ After this paper was submitted, Hailey & Mori (2002) proposed that the features 
presence of substructures in the Unes. 



require a very small initial kinetic energy, since at this distance 
the SNR would be ^700 pc above the galactic plane, where the 
intersteUar density is expected to be very small. In conclusion, 
we believe that a distance of ~ 2 kpc and an age of the order of 
^ 10"* years are more likely values. 

Although we cannot rule out other possibifities (e.g. one or 
more glitches in the last two years, significant timing noise, or a 
faint binary companion), the simplest explanation to reconcile 
Tc with a young age is that IE 1207-5209 was born spinning 
at a period close to its current value. 

Thus IE 1207-5209 appears as a unique example of a young 
pulsar with a small rotational energy loss. Its rotational energy 
loss, Eroi^ 10^"* erg s"', is the smallest one of all the radio/X- 
ray pulsars with reliable SNR's associations (Kaspi & Helfand 
2002). Such "Crab-like" or "Vela-like" pulsars have typically 
Erot > 10^^ erg s~' and X-ray luminosity in the range 10^^-10^' 
erg s~' (see, e.g., Possenti et al. 2002). Indeed, no plerion has 
been detected either in radio or X-rays, nor any hint has been 
found of 7-ray emission from our object. On this ground, the 
comparison with the similar age Vela pulsar is striking: beside 
being active as a radio pulsar, Vela is also the brightest source 
in the sky at E>100 MeV, and is powering a bright radio/X- 
ray synchrotron nebula. If the lines are the signature of the 
presence of metals in the atmosphere, their absence in the older 
objects mentioned above can be explained by a small layer of 
hydrogen accreted, even at a modest rate from the interstellar 
medium over a time span of a few hundred kyrs. The absence 
of lines (or the difficulty to detect them) in young but energetic 
pulsars fike Vela could be due to the possible perturbing effects 
of the non-thermal particles accelerated in the magnetosphere. 



4. CONCLUSIONS 

The XMM-Newton observations reported here show unam- 
biguously that the absorption features independently discovered 
with Chandra (Sanwal et al. 2002) in the X-ray spectrum of 
IE 1207-5209 are phase-dependent and have significant sub- 
structure. This supports an interpretation in terms of atomic 
transitions in regions with different temperature and magnetic 
field on the neutron star surface. A detailed analysis of phase- 
resolved spectra with higher statistical quality will undoubtedly 
provide important information on the surface composition and 
other neutron star parameters. 

The period value reported here is identical, within the uncer- 
tainties, to that obtained two weeks later (Pavlov et al. 2002a), 
but our smaller error reduces the uncertainty on the pulsar spin- 
down rate P=(1.98±0.83)x 10"'"* s s"*. Since the association 
with the relatively young SNR remnant is extremely Ukely, we 
consider the high characteristic age of IE 1207-5209 as evi- 
dence that this pulsar was born with a long period. 

We suggest that the combination of a young age and a low 
level of magnetospheric activity are favorable conditions for the 
detection of atomic spectral features from Z> 1 elements, which 
would be either blanketed by a thin layer of accreted hydrogen 
in older objects or masked by non-thermal processes in young 
energetic pulsars. If this scenario is correct, the most promising 
candidates for fine detection are the central X-ray sources in 
young SNRs which do not show synchrotron nebulae, like Cas 
A (Mereghetti, Tiengo & Israel 2002), Puppis A (ZavUn, Trum- 
per & Pavlov 1999), and G 266.2-1.2 (Pavlov et al. 2001b). 

in IE 1207-5209 are produced by He-like Oxigen or Neon and predicted the 
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Table 1 

Results of PN fits including absorption lines 



Blackbody Blackbody + Power law H atmosphere 
+ 2 Gaussians + 2 Gaussians + 2 Gaussians 



(10^° cm-2) 


3-0'^i 2 


10+4 


9+2 


kT (keV) 


0.244±0.005 


0.219±0.008 


0.096±0.005 


Rfg (km) 


I.7±0.2 


2.1±0.3 


13±1 


Olph 




2.9!?4 




El (keV) 


0.74±0.02 


0.73^«i 


7 1+0 02 


cri (keV) 


0.13±0.03 


0.12±0.04 


0.12±0.04 


EWx (eV) 




Qo+30 
^^-120 


_Q1+15 


E2 (keV) 


I.36±0.03 


1.36±0.02 


1.37±0.02 


CT2 (keV) 


0.10±0.05 


0.08±0.04 


0.07±0.03 


EWi (eV) 


-107^36 


-671^ 


-6318 


Fol3fev(ergcm-2s-i) 


1.96x10-12 


1.96x10-12 


1.95x10-12 


Fp^ (erg cm"-^ s~') 




5.5x10-" 




Lf (ergs-') 


1.15x10^3 


1.17x10^^ 


1.8x10^3 


x'/dof 


1.54 


1.10 


1.16 


dof 


75 


73 


75 



" Radius at infinity for an assumed distance of 2 kpc. 
* Observed flux. 

Observed flux of the power law component only, 0.3-3 keV. 

Bolometric luminosity (excluding the power law component) for d=2 kpc. 
All the errors are at the 90% c.l. for a single interesting parameter 
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Table 2 

Results of phase resolved spectroscopy 





A 




B 


C 


D 




0.8-1. 




0-0.25 


0.25-0.55 


0.55-0.8 


Nh (lO^" cm-2) fixed 


3.0 




3.0 


3.0 


3.0 


kT (keV) 


0.241±0.004 


0.250±0.004 


0.245±0.004 


0.242±0.004 


41 (km) 


1.8 




1.7 


1.6 


1.7 


El (keV) 


0.78±0.04 


0, 


.76±0.04 


0.74±0.02 


0.71 ±0.02 


CTi (keV) 


0.12±0.02 


0, 


.16±0.04 


0.13±0.02 


0.08±0.02 


rf 


0.29±0.08 


0, 


.28±0.08 


0.35±0.07 


0.37±0.11 


El (keV) 


1.36±0.04 


1, 


.37±0.03 


1.37±0.02 


1.35±0.02 


(72 (keV) 




0, 


.18±0.04 


0.09±0.02 


0.06±0.02 


r2 


0.19±0.09 


0, 


.30±0.08 


0.52±0.13 


0.56±0.16 


X^/dof 


1.39 




1.00 


1.48 


1.31 


dof 


63 




68 


69 


67 



" Radius at infinity for an assumed distance of 2 kpc. 

* Relative line depth = 1. - F (E;,>,e) / Fc (E/,„e), where F is the total flux and Fc is the flux of the continuum only 
AH the errors are at the 90% confidence level for a single interesting parameter 



